Two regimes for effects of surface disorder on the zero-bias conductance peak of 
tunnel junctions involving rf-wave superconductors 
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Impurity-induced quasiparticle bound states on a pair-breaking surface of a d-wave superconduc- 
tor are theoretically described, taking into account hybridization of impurity- and surface-induced 
Andreev states. Further a theory for effects of surface disorder (of thin impurity surface layer) on 
the low-bias conductance of tunnel junctions is developed. We find a threshold Uc for surface impu- 
rity concentration ns, which separates the two regimes for surface impurity effects on the zero-bias 
conductance peak (ZBCP). Below the threshold, surface impurities do not broaden the ZBCP, but 
effectively reduce its weight and generate impurity bands. For low ns impurity bands can be, in 
principle, resolved experimentally, being centered at energies of bound states induced by an isolated 
impurity on the surface. For larger ns impurity bands are distorted, move to lower energies and, be- 
ginning with the threshold concentration ns ~ nc, become centered at zero energy. With increasing 
ns above the threshold, the ZBCP is quickly destroyed in the case of strong scatterers, while it is 
gradually suppressed and broaden in the presence of weak impurity potentials. More realistic cases, 
taking into account additional broadening, not related to the surface disorder, are also considered. 

PACS numbers: 74.45.-l-c, 74.72.-h 
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Inhomogcncitics in anisotropically paired supcrfluids 
and superconductors can noticeably modify properties of 
adjacent regions on the scale of the coherence length. In 
particular, quasiparticle Andreev bound states can arise 
in those regions near impurities, surfaces or interfaces. 
Zero-energy Andreev bound states near impenetrable 
smooth surfaces of d-wavc superconductors arc known to 
be the signature of a sign change of the order parameter 
in a quasiparticle reflection event. In high-temperature 
superconductors the zero-energy states are identified, in 
particular, as resulting at low temperatures in the zero- 
bias conductance peak (ZBCP) of NIS-junctions (see, for 
example, review articlesi*^ and references therein) . 

In the present paper we report a theory for effects of 
surface disorder on the ZBCP and, more generally, the 
subgap conductance. Our starting point is the problem 
of quasiparticle bound states induced by an isolated im- 
purity, situated close or directly on (110) surface of a 
d-wave superconductor. Spectra of the bound states can 
be found from the equation for the poles of the t-matrix, 
which takes the form 
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Here ri,„p is the impurity position and Vimjy the weight of 
the potential V{r) = Vimp5{r — Timp)- For simplicity, we 
disregard here the effect of the impurity on the order pa- 
rameter. The Nambu matrix retarded Green's function 
Go{r, r'; e) describes the system in the absence of the im- 
purity. The quantity Go{rimp, Vimp] e) with identical first 
and second coordinates is directly associated with the 
quasiclassical Green's function averaged over the Fermi 



surface: Go(r". 
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Here we omitted the term, which should be eventually 



included into the renormalized chemical potential; Nf 
is the normal state density of states on the Fermi sur- 
face per one spin direction, is the respective Pauli ma- 
trix in particle-hole space. In the presence of (110) or 
(100) surfaces in a c?-wave superconductor, the averaged 
off-diagonal components of the Green's function vanish 
{foiPf,x;e))^^ = (/o+(py,x;e))^^ = throughout the 

superconductor and for any e. In this case we obtain the 
following equation for the spectrum of impurity bound 
states: 
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where dimensionless impurity potential is introduced u = 
T^VimpNj. Near (100) surface the quasiclassical Green's 
fimction of a d2.2_j,2-wave superconductor does not de- 
pend on spatial coordinates and takes the same form as 

in the bulk: gQ{pf,r]e) = —ine j — |A(pj)p . Av- 
eraging this expression over a cylindrical Fermi surface 
and substituting into Eq. Q leads to well established 
results on energies of impurity bound states in d-wave 
superconductors^. Qualitatively, a well pronounced 
low-energy resonance forms in the case of sufficiently 
large strength of the impurity potential, for instance, 
close to the unitarity limit. No resonance is present in 
the case of weakly scattering impurity, when subgap im- 
purity states lie mostly at higher energies, being strongly 
broaden. 

We demonstrate below, that impurity states near (110) 
surface have quite different structure as compared with 
the (100) case, as an effect of a hybridization of impurity- 
and surface-induced Andreev states. In the presence 
of (110) surface the quasiclassical Green's function of a 
(ij.2_y2-wave superconductor is spatially dependent and 
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has a pole at zero energy, associated with zero-energy sur- 
face Andreev bound states. For this reason the spectrum 
of impurity bound states becomes dependent on the dis- 
tance of the impurity from the surface. Assume, for sim- 
pUcity, that the impurity is situated quite close or directly 
on (110) surface, so that the zero-energy surface states 
have strong influence on the formation of the impurity 
states. In contrast with the (100) case, a well pronounced 
low-energy impurity resonance Eimp = e'imp + i^'lmp^ 
k'j'npl ^ l^'impl takes place in this case only for weakly 
scattering impurity potential. At low energies the pole- 
like term dominates the Green's function and in the 
first approximation ((;o(p/, a; = 0; e))^^ = j-pjg, where 

r = 7r(|A(pj:)|)^^ is the residue of the pole-like term 

averaged over the Fermi surface. The quantity |A(p/)| 
coincides with the modulus of the bulk order parame- 
ter, if one disregards the surface pair-breaking. In gen- 
eral, it is a normalized integral over the superconductor, 
containing a spatially dependent (self-consistent) order 
parametep^. Inserting the pole-like term into the equa- 
tion for impurity bound states, we find the energy of 
the resonance for weakly scattering impurity potential 
{u <^ 1): e^mp = ±ur /-K . In order to find the decay 
rate, one should keep also the first low-energy term in the 
imaginary part of the Green's function, averaged over the 
Fermi surface. For nonzero low energies, narrow vicini- 
ties near the nodes of the order parameter dominate the 
averaging of the imaginary part. Assuming a cylindrical 
Fermi surface, we find {go{pf,x ~ 0; e))^^ ~ eTzO 
where Ai is the slope of the order parameter at the node, 
taken in the bulk: Ai = \A'{(pnode)\- With this low- 
energy expression for the averaged Green's function we 
obtain from Eq.lQ for u ^ 1: 
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±ur/'K - ir^u-^/in^ Ai) . (3) 



Self-consistent calculations give for (110) surface orien- 
tation r « Aq with a good accuracy, where Aq is the 
maximum value of the bulk order parameter A(py) = 
Ao sin 2(p. Impurity-induced subgap states in d-wavc su- 
perconductors have Andreev origin. An imaginary part 
of the energy e'-„^p is associated with an escape of the 
quasiparticle from the impurity into the bulk of the su- 
perconductor. This can take place due to impurity scat- 
tering, if for an acquired momentum direction the quasi- 
particle energy exceeds the modulus of the c?-wave order 
parameter. For sufficiently small energy |e^„jp| <C Aq 
a quasiparticle can run away from the impurity only 
with momenta in a narrow vicinity of nodal directions of 
the order parameter. Then a well pronounced impurity- 
induced quasiparticle resonance arises. For a larger sub- 
gap energy of impurity-induced state, the escape can oc- 
cur in a wider region of momentum directions and the 
resonance is more broaden. As seen from Eq.J^)), the 
resonance is well defined for weakly scattering impurity 
potential u <C 1. In the opposite case of large strength 
of impurity potential 1, the impurity states formally 



Similar to the impurity in the bulk of d-wave 
superconductor—, the local density of states (LDOS) 
at the impurity site near (110) surface manifests the 
particle-hole asymmetry: for repulsive impurity poten- 
tial a resonance peak in the LDOS at this site takes 
place only for a; < 0, while for attractive potential at 
a; > 0. This property is directly associated with van- 
ishing anomalous Gor'kov Green's function, describing 
the d-wave superconductor with (110) surface in the ab- 
sence of the impurity and taken at one and the same site 
Fo{r,mp,rimp,£) = Nf (/o(p/, n^^; e))^^, = 0. 

The resonance energy is a function on the distance 
Ximp between the impurity and the (110) surface. At 
sufficiently large distances Ximp 2> i ^^'^o and not neces- 



sarily for weak scatterers, we obtain 
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The weight of the corresponding 
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peak in the LDOS at the impurity location diminishes 
cx x~^p with increasing Ximp, similar to the weight of 
the zero-energy peak originated from the surface bound 
states. 

It turns out that changes in the LDOS due to the im- 
purity resonance, take place entirely at the expense of 
the zero-energy surface states. The presence of a single 
impurity near (110) surface induces four quasiparticle im- 
purity states: two spin degenerated states with positive 
and two with negative energies. Their contribution to 
the integral DOS is 2, since each spin-polarized Andreev 
state with nonzero energy contribute one half to the inte- 
gral DOSiS. In its turn, the number of the zero-energy 
Andreev surface states in the presence of one impurity 
reduces exactly by two (by one per spin polarization), 
as compared with the case of no surface defects. Thus, 
the total number of the surface and the impurity quasi- 
particle states keeps constant. This is valid, actually, for 
any number of impurities at (or near) the surface. The 
resonance energy and decay rate depend also on surface- 
to-crystal orientation. The above consideration can be 
easily generalized, with minor modifications, to a wide 
range of the orientations. Only for surface orientations 
very close to (100), when the residue of the zero-energy 
pole-like term becomes negligible (due to a small frac- 
tion of trajectories where the zero-energy surface states 

take place), a small low-energy term — ^ In dom- 
inates the averaged Green's function. As a result, the 
spectrum of impurity bound states on (100) surface coin- 
cide with the case of impurity-induced states in the bulk 



of a d 



-wave superconductor. 
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Aq and the resonance itself is ill de- 



Consider further the low-energy LDOS and the low- 
bias tunnel conductance for the d-wave superconductor 
with (110) surface in the presence of many impurities. 
Since the conductance is sensitive mostly to impurities 
in the vicinity of the surface, we introduce below a thin 
surface layer of strongly disordered normal metal, which 
differs from the adjacent only by its very short mean 
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free path I. We consider a clean d-wave superconductor 
in the half-space x > d. The surface layer with impu- 
rity concentration Uimp is placed at < x < c? between 
a smooth impenetrable specularly reflecting surface at 
X ~ and a clean superconductor at x > d. There 
are no additional potential barriers at a; = d, between 
the disordered layer and the superconductor. The thick- 
ness d is assumed much less than all other characteristic 
scales in the problem. The model is very similar to the 
thin dirty layer model (TDL)S except that it incorporates 
the self-consistent t-matrix approach, allowing for arbi- 
trary strength of quasiparticle scattering by impurities, 
while the TDL considers only the Born scattering. The 
other difference is that real surface layer of normal metal 
with isotropic impurities is introduced here, while special 
anisotropic scatterers (Ovchinnikov impurities) are used 
in the models of surface roughness^iifl. Our results, pre- 
sented below, can be easily reformulated for the case of 
Ovchinnikov impurities as well. 

The quasiclassical 2x2 particle-hole matrix retarded 

9 I 

quasiparticle excitations and obeys Eilcnberger's equa- 
tions. For thin disordered surface layer with the impurity 
self-energy aimpix;£) ~ (Timp(0;e), the surface impurity 
concentration ns = riimpd and the surface impurity self- 
energy = 2iaimpix = 0;e)d enter the final results. 



Green's function g{pf,x;e) = 



describes 



vanish 



Off-diagonal elements of ({e) and (^g{pj:,x 

in the case of (110) surface, so that C(e) = C{e)T3. 

The shape of the zero-energy peak in the density of 
states depends, to a great extent, on the surface self- 
energy C{e). The surface self-energy, in its turn, is con- 
trolled by the Green's function averaged over the Fermi 
surface. We admit isotropic scatterers with arbitrary 
strength of potentials. Introducing a dimensionless sur- 
face impurity concentration ns = 2ns/{TThvfNf), the 
surface self-energy parametrizes as: 



ivfus (.9(P/>0;e))g^ 



(4) 



For the low-bias conductance only the Green's function 
on the surface at low energies is of importance. The 
analytical solution we found for the model gives for low 
energies 



g{pf, X = 0; e) — —in coth 



C(e) 



v/|cos(p| |A(p/)| 



(5) 

Eqs. I@J, lO were derived exploiting the self-consistent 
spatial profile of the order parameter. Our further goal 
is to solve them jointly in some important particular cases 
and describe analytically the respective LDOS as a func- 
tion of the energy, the surface impurity concentration 
and the strength of impurity potential. Also we repre- 
sent below our fully self-consistent numerical results for 
conductance spectra, based on quasiclassical equations 



and boundary conditions for Riccati amplitudes of the 
quasiclassical matrix Green's functioniiii^. 

It turns out that the impurity layer is unable to 
broaden zero-energy surface bound states, if surface im- 
purity concentration is less than the characteristic thresh- 
old: C(0) = for fis < fic ~ {\cosLp\)g^. As we 
found analytically, the delta-peak at zero energy occurs 
in the LDOS for all surface impurity concentrations be- 
low the threshold ns < Uc (see, e.g., Eq. O). Al- 
though, the weight of the peak is quite sensitive to ns 
(see below). For the simplest model of a quasi- two- 
dimensional superconductor with cylindrical Fermi sur- 
face he = 2/7r. The value he is determined by the sym- 
metry of pairing and does not depend on the particular 
spatial profile of the order parameter and the form of 
its basis functions. The threshold concentration of sur- 
face impurities exactly coincides with the number A^o(O) 
of the zero-energy bound states per unit area of (110) 
surface in the absence of the disordered surface layer: 
Nq{0) = Uc = {TT/2)NfV fhc- The effect is associated 
with detaching the states from the zero-energy peak in 
forming impurity bands. To a certain extent, this bears 
an analogy to what is known for impurity broadening of 
Landau levels and lifting of their degeneracy in the two- 
dimensional electron gas in a strong magnetic fieldi^ii^*i^. 
In the presence of the disordered layer the number of 
zero-energy surface states linearly depends on the surface 
impurity concentration hs and vanishes at the threshold: 
A^o("-s) = No{0){hc — hs)/hc- The threshold concen- 
tration remains equal to the number of the zero-energy 
states on the clean surface also for arbitrary surface mis- 
orientation angle a, not only for a = 45°. For cylindrical 
Fermi surface nda) = NfVfy^l — |cos2a|. It approaches 
its maximal value for a = 45°. 

The behavior of the LDOS qualitatively differs below 
and above the impurity concentration threshold. In the 
presence of the surface impurity layer impurity bound 
states transform into two impurity bands, one with pos- 
itive and the other with negative energies. For low con- 
centrations of the Born impurities the bands are situated 
at low energies and the LDOS for the impurity bands 
takes the form 



-is) = ^v'(W^/2)^-(e-4„,)^ 



(6) 



This expression is valid in the energy range, where the 
argument of the square root function is positive. The 
characteristic bandwidth W ^ AoU\/ns and the height 
of the LDOS in the center of the band Vmax ~ Vns/u 
arc introduced in Eq. Numerical coefhcients in ex- 
pressions for W and v depend only on |A(p/)|. We 
omit here respective cumbersome analytical formulas. In 
a good agreement with self-consistent calculations, the 
center of the Born impurity band is nothing but the po- 
sition of the impurity resonance on the (110) surface: 
\e'^jyjp\ = u(|A(p/)|)g^ . As seen from these estimations, 

the impurity hump rises, becomes more narrow and goes 
towards the zero energy, when the potential strength de- 




FIG. 1: Tunnel density of states on (110) surface as a func- 
tion of energy. Solid line represents tunnel density of states 
for a smooth clean surface. Other curves describe effects of 
the surface impurity layer with various strengths of impu- 
rity potentials and surface impurity concentrations. Intrinsic 
broadening 0.002rc is introduced for resolving the zero-energy 
J-peak. 



FIG. 2: Tunnel density of states on (110) surface as a function 
of energy, taken for various surface impurity concentrations 
below and at the threshold. The strength of the impurity 
potential is identical for all curves: u — 1. Inset shows the 
tunnel density of states for impurity concentrations above the 
threshold: solid line - ns = 2.5nc, u = 0.3, dashed line - 
ns = 1.2nc, M = 1, dotted line - ns = 1.07nc, u = 10. 



creases. Eq. describing the shape of the impurity 
band in the LDOS, is vahd only for small surface impu- 
rity concentration Wmax ^ \/ns ^ 1 and implies that 
the width of the impurity band significantly exceeds the 
width of separate impurity resonance W ^ l^impl- 

Fig. ^ displays the energy dependence of the tunnel 
density of states (i.e. the conductance at T = 0) for clean 
(110) surface, as well as for surface disordered layers with 
various strengths of impurity potentials. Here and below 
the momentum dependence of the transparency of NIS 
junction is taken in the simple form D[tf) ~ Dq cos^ Lp, 
suitable for thin and high potential barriers. Tunnel 
DOS and conductance are normalized to their normal 
state values. At those energies, where the impurity peak 
stands out well against the background, the tunnel den- 
sity of states is simply proportional to the LDOS on the 
surface with the energy-independent coefficient of the or- 
der of unity. Eq.® describes with a good accuracy im- 
purity peaks on two curves with u = 0.3 in Fig.^ The 
impurity peaks are situated sufficiently close to the zero- 
energy peak. In the case u = 0.3, ns — 0.2nc the impu- 
rity peak has a small asymmetric distortion, which is not 
taken into account in Eq.JBJ. With increasing u the im- 
purity bands become less pronounced, more asymmetric 
and broaden. For sufficiently large u impurity bands are 
strongly broaden at low ns over all subgap region. 

Fig.[21shows the evolution of impurity bands with vary- 
ing the surface impurity concentration below the thresh- 
old for the impurity potential u = I. The impurity 



peak goes towards the lower energies and rises when the 
concentration increases. Beginning with the threshold 
concentration impurity bands are centered on zero en- 
ergy, for any strength of impurity potential. At the 
threshold concentration an impurity-induced zero-energy 
peak in the tunnel DOS arises with a specific power-law 
spectral shape vd{£) « l^l"^^'^ in the low-energy region 
ln(Ao/|e|) ^ 1. This peak takes place for arbitrary 
strength of impurity potential. 

The weights of the ZBCP and the zero-energy peak 
in the LDOS, associated with the surface zero-energy 
Andreev states, reduce with increasing ns- Below the 
threshold they vary linear in [he — ns), within the log- 
arithmic accuracy, and vanish at ns = nc- Close to the 
threshold < fic — ns <^ nc the low-energy momentum 
resolved LDOS 



7r| cos ip\ {fie - ns) 



So In 
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O^max 

ns 



(7) 



Here 5(.t) = (cos^ <p S{x-\A{pf)\)) Sq = S{x +0), 

5oln5 = S' (x)\vL{K,naxlx)dx and A^aa is a maxi- 
mal value of |A(p/)| over Fermi surface. For cylindrical 
Fermi surface Sq = 2/7rAi. With spatially constant or- 
der parameter we also find iS = 2. Eq. (jZJ is valid for low 
energies |e| <^ {he — hs) Ai, {he — hs)^^'^uAi, where one 
can disregard the contribution from impurity bands. 
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Above the threshold hs > fie the solution of Eqs.Q), 
results in the following angular resolved surface den- 
sity of states at zero energy 



3^ 

2m2 



In 



13.1 



ns ■ 



1/2 



^6(ns - fic) 



I cosv^l 



(8) 



This finite zero-energy value of the density of states, to- 
gether with the low-energy corrections (see, e.g., Eq.®), 
describes the broaden zero-energy peak above the thresh- 
old. The condition C(0)/vj <C 1 implied in the derivation 
of Eqs.® applies not only close to the threshold. For suf- 
ficiently weak scatterers, when the large term ZtP' jlii^ 
dominates the logarithmic function, the expression for 
C(0)/v/ remains also small for hg ^ ric- The Born ap- 
proximation works in this case and we find v{j^^ x = 
0;e = 0) = I cos</3|/(pric)"^/^, where p = v?hs = d/l. 
Here I = VfT and r are the mean free path and the re- 
laxation time in the disordered layer. The inset of Fig.[2| 
demonstrates the behavior of the low-energy tunnel den- 
sity of states above the threshold. The heights of ZBCP, 
which can be obtained from Eq. ©, are in agreement 
with those in the inset of Fig. [3 within few percents. 

Usually, the Born approximation is justified, if the first 
term in the denominator of Eq. dominates the second 
one, i.e. the relation li^j (^(p/, 0; e))^^ I/tt^ <C 1 satisfies. 
The unitary approximation works in the opposite limit. 
In ordinary normal metals | ((7(p/, 0; e))^^ I/tt^ = 1 and 
the condition for the Born (unitary) approximation to be 
valid reduces to the standard form: <C 1 {u^ ^ 1). 
In superconductors the quantity | {g{pf,0;e))g^ \/n^ can 
strongly deviate from unity at low energies. In the ab- 
sence of zero-energy (low-energy) states we have at low 
energies | (g(p/, 0; e))^^ I/tt^ ^ 1. This expands the re- 
gion of applicability of the Born approximation in su- 
perconductors at low energies. However, in the pres- 
ence of the zero-energy surface bound states the pole- 
like term dominates the Green's function and takes very 
large values on the surface at low energies. Moreover, 
since the zero-energy states are dispersionless states, the 
Green's function, taken on the boundary at low ener- 
gies, remains to be quite large even after averaging over 
the Fermi surface. Then, at sufficiently low energies we 
obtain | (p(p/, 0; e))^^ I/tt^ :$> 1. This condition strongly 

modifies the behavior of the results, obtained in the Born 
and in the unitary approximationsi^. It is also seen, that 
even though there are Born impurities in the disordered 
surface layer on the normal metal sample {u^ ^ 1), in 
the case of the superconducting substrate the presence 
of zero-energy surface states places additional restric- 
tions on admitting the Born approximation to the same 
scatterers. They can scatter low-energy quasiparticles 
near the surface with a large effective strength. This 
takes place, if the impurity concentration is not suffi- 
ciently high and the zero-energy peak in the LDOS is well 
pronounced iiZi Under certain conditions the weak impu- 



rity potentials <C 1 can be described even with unitary 
approximation in the presence of the zero-energy surface 
states. This applies below the threshold to describing 
the zero-energy peak and turns out to be very important 
for our results. This is very different from various other 
circumstances where the Green's function doesn't take 
large values and the Born approximation can always be 
justified for sufficiently weak scatterers. 

We find, that the Born approximation applies for de- 
scribing sufficiently high zero-energy peaks in the den- 
sity of states, only if the surface impurity concentra- 
tion well exceeds the threshold surface concentration: 
ns 3> JT-c^ This condition can be quite restrictive. A 
simple estimation gives Uc ^ hvfNf an^, a being the 
interatomic distance and rig the electron concentration. 
For optimally doped YBa2Cu307_5 cut with (110) sur- 
face, we evaluate the threshold surface concentration as 
ric « 1.3 X 10^'*cm~2. This is sufficiently large (although 
well accessible) value of the threshold for the surface con- 
centration of point defects. 

Conditions for the applicability of the Born and the 
unitary approximations modify further for surface im- 
purity concentrations close to the threshold. The self- 
consistent analysis of Eqs. Q), jSJ shows that at ns « ric 
the two large terms, which usually dominate the low en- 
ergy behavior of (g{pf, 0; i almost cancel each other. 
Then conventional criteria for the Born (m ^ 1) and the 
unitary (u 3> 1) approximations recover with a logarith- 
mic accuracy. In particular, this kind of conditions for 
the Born or the unitary approximations arises in Eq.®, 
depending on what term in the numerator dominates. 

A characteristic feature of the zero-energy peak in 
the LDOS above the surface concentration threshold, is 
a cusp-like shape of the peak. Indeed, the first low- 
energy correction to the LDOS turns out to have a 
characteristic energy dependence ((5^/(^/^, a: = 0; e))^^ = 
(1^(93, 0, e) — z^((p, 0, 0))_5^ oc |e|. The cusp-like shape of 
the broaden zero-energy peak is clearly seen on the in- 
set of Fig. 13 (see also Fig. (6c) in Ref. [ij, obtained with 
the thin dirty layer model and Ovchinnikov impurities). 
While the zero-energy peak in the LDOS originates from 
the sign change of the order parameter in a quasiparticle 
reflection from the surface, the cusp-like behavior of the 
peak comes from the narrow vicinities of the order param- 
eter nodes. Since for (110) surface only the momentum 
directions close to the surface normal are of importance, 
the tunnel DOS at low energies acquires the same cusp- 
like behavior ~ |e|, although with a different numerical 
coefficient. In the unitary limit, when the zero-energy 
density of states is described by Eq.© with u 1, we 
get 



((5z^(^,0;e))_5^ 



^'(0) 



8 i/2(o)(^5_^^)_ 4/(37^3)- 



This term becomes comparable with (z^((^, x ~ 0; 0))^^ 
1^(0) for e ~ (ns — n^^"^ 



(9) 



Our estimations for suffi- 



ciently high and narrow ZBCP show that the width of the 
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peak manifests analogous dependence 7 ~ {fis — fic)'^^^ ■ 
In the Born limit we obtain {Siy(ip,x = 0;e)) g^, ~ 
--7r|£|/(4pAi). This quantity is of the order of the 
zero-energy surface density of states ~ l/lpficY^^ for 
e '--^ 7 cx p^/^. For the Born surface disordered layer 
({0)/vf = \/ pfic- It follows from here 7 (x ^ oc r~^/^. 
We note, that the relation 7 cx r~^/^ is a common fea- 
ture of the broadening of the zero-energy surface bound 
states by Born scatterers, doesn't matter whether they 
are spread out as bulk impuritiesi^, occupy a surface 
layer with the thickness of the order or larger than the 
coherence length^S, or are collected only in a very thin 
surface layer as it is the case in the present paper. The 
width of a well pronounced impurity band turns out to 
be (X y^riimp also below the threshold, as it is for W in 
Eq.®, as well as in some other circumstancesSLSS. 
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FIG. 3: LDOS on (110) surface as a function of energy, taken 
for three values of strength of impurity potential u and for 
various bulk impurity concentrations. Dimensionless concen- 
tration of bulk impurities hf, is defined as tij, = 2Tr'^NfTcnt,. 
The lower curve corresponds to clean superconductor ht = 0. 
With hb = 0; 0.001; 0.002; 0.005; 0.01; 0.02 one goes from 
the lower to the upper curve on each panel. Impurity humps 
are present for u — 0.3, u — 1.0 for a certain range of concen- 
trations rit- 

Consider further more realistic cases, taking into ac- 
count also some other reasons, not related to the surface 
disorder, which can result in a broadening of the ZBCP in 
both regimes studied above. Let, for example, impurities 
occupy not only the surface layer with surface impurity 
concentration ns, but also are spread out over an exten- 
sive region of a superconductor as bulk impurities with 
small bulk concentration n^. Since energies of surface im- 
purity states depend on distance Ximp, in the presence of 
bulk impurities it is not obvious whether impurity bands 
are fully smeared out and merge in a periphery of the 
zero-energy peak, or they still survive, being centered at 
energies of surface impurity states Eq.(|3Jl for low ns- For 
making clear this question we have studied, firstly, the 
effect of bulk impurities only. Fig. O shows the LDOS 
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FIG. 4: LDOS on (110) surface in the presence of bulk and 
surface impurities, taken for various bulk impurity concentra- 
tions nb, as a function of energy. The strength of impurity po- 
tential It = 1, the surface impurity concentration ns ~ 0.2nc. 

on (110) surface for a superconducting d-wave half-space 
with small bulk impurity concentration rif,, which is con- 
stant everywhere up to (110) surface. As seen, surface 
impurity bands survive in the LDOS for weak scatterers 
(u = 0.3 and u = 1) in a range of bulk impurity concen- 
trations rif,. Although, the bands arc less pronounced and 
take place in a more narrow region of ri;,, as compared 
with Figs.m^lfor the surface impurity layer model. With 
increasing potential strength or/and impurity concentra- 
tion 71b, impurity bands are gradually smeared out and 
merge in the zero-energy peak. 

Effects of bulk impurities in the presence of thin sur- 
face impurity layer on the LDOS on (110) surface, are 
shown in Fig. 0] for ns — 0.2nc and various Ub- Di- 
mensionless concentration of bulk impurities is defined 
as rib = 2TT'^NfTchb. One can see, that bulk impuri- 
ties control the broadening of the zero-energy peak, when 
surface disorder lies well below the threshold. The zero- 
energy peak and the impurity bands fully merge into one 
broaden peak for comparatively small concentration of 
bulk impurities Ub ~ 0.047r^iV/rc. Even in this case the 
height of the zero-energy peak is sensitive to the sur- 
face disorder, which controls relative weights of central 
and peripheral regions of the broaden zero-energy peak. 
The "fine structure" of the peak can be resolved only 
for less rib- Then well above the threshold concentration 
the surface layer controls the broadening. Fig. [3 shows 
the evolution of the low-energy LDOS with the surface 
impurity concentration ns in the surface layer, in the 
presence of bulk impurities {nb = O.ObnNfTc). The inset 
of Fig. [31 shows the broadening 7 of the zero-energy peak 
as a function of surface impurity concentration ns/nc, 
taken for various bulk impurity concentrations hb. The 
broadening 7 is the half of the width of the zero-energy 
peak: 1/(0) — 22/(7). As seen in the inset of Fig. |S1 al- 
ready at small concentrations bulk impurities transform 



7 




FIG. 5: LDOS on (110) surface in the presence of bulk and 
surface impurities, taken for various surface impurity con- 
centrations ns- The strength of impurity potential u = 1, 
the bulk impurity concentration hb = 0.025. Inset shows the 
broadening 7 of the zero-energy peak as a function of surface 
impurity concentration ng/nc, taken for various bulk impu- 
rity concentrations hi,. 



an abrupt change of the regimes of the broadening at 
the threshold into a gradual crossover from one regime 
to another. 
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FIG. 7: Tunnel conductance on (110) oriented surface as a 
function of applied voltage for different impurity concentra- 
tion ns = 0.2nc, ns = 0.4nc, ns = 0.7nc, ns ~ l.Oric, 
ns = lAnc- Higher ZBCP corresponds to lower concentra- 
tion. Temperature and scattering impurity potential are the 
same for all curves; T — 0.04rc, u = 1. Inset shows the 
broadening 7 of the ZBCP as a function of surface impurity 
concentration, taken for various temperatures: solid line - 
T = 0, dashed line - T = OMTc, dotted line -T = OMTc. 
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FIG. 6: Tunnel conductance on (110) surface for various tem- 
peratures as a function of the applied voltage. With increas- 
ing the temperature the ZBCP and impurity states broaden 
and strongly overlap. The strength of impurity potential 
and the concentration are identical for all curves: m = 1, 
ns ~ 0.2nc. 

Effects of finite temperatures also contribute to a 



broadening and a suppression of the ZBCP, does not mat- 
ter in the presence or absence of surface disorder. The 
tunnel conductance as a function of the applied voltage 
is shovifn on Fig. Elfor various temperatures. The tem- 
perature broadening of both the zero-energy and the im- 
purity peaks takes place. The peaks fully merge into the 
one broaden peak above a comparatively low tempera- 
ture, which is ~ 0.06Tc for the particular case, consid- 
ered on Fig. El The "fine structure" of the ZBCP arises 
below this temperature. In the cases we consider, impu- 
rity bands result in two humps (small satellites) whose 
heights are usually noticeably less than the ZBCP in the 
self-consistent calculations. Only for concentrations very 
close to the threshold the impurity peaks can become 
high and centered very close to zero energy, while the 
zero-energy peak is already exhausted. We tried to single 
out, whether this kind of "spontaneous splitting" of one 
broaden ZBCP into two separate low-energy peaks can be 
observed below a characteristic temperature in the case 
of high spectral resolution. Within the self-consistent t- 
matrix approximation and for all sets of parameters we 
used, three peaks always merge into one broaden peak 
before a "spontaneous splitting" described above could 
show up. 

Fig. [7| shows the temperature broadening of the ZBCP 
for various surface impurity concentrations. As seen in 
the range of low ns, the broadening of the ZBCP, being 
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fixed by the temperature, is almost insensitive to ng. At 
the same time the height of the peak noticeably decreases 
with increasing ng. The inset of Fig.|Z|demonstrates that 
effects of finite temperatures, analogously to the effect 
of bulk impurities, transform the abrupt change of the 
regimes at the threshold into a gradual crossover from one 
regime to another. The surface impurity concentration 
can strongly influence the width of the ZBCP only above 
the crossover region, while the height of the ZBCP can 
be sensitive to the surface disorder for any 715. 

Results of the present paper apply to tunnel junc- 
tions, when the quasiparticle escape from the supercon- 
ductor to the normal metal does not have strong influ- 
ence on the broadening of the zero-energy state o^i^? . For 
higher transparencies this mechanism of the broadening 
becomes important. An information to what extent im- 
purities and effects of finite transparency are involved 
into forming the broadening of the ZBCP in a particular 
NIS junction, in principle, could be obtained experimen- 
tally, e. g. by measurements of the differential shot noise 
at low voltages^^. 

Thus, in the case of strong broadening of the ZBCP, 
induced by sources, which are not related to the sur- 
face disorder, impurity bands and the zero-energy peak, 
associated with Andreev surface states, merge into one 
broaden ZBCP. This can be associated also with the ex- 
perimental resolution of the spectral weight. Even in 
this case the height of the ZBCP is quite sensitive to 
the surface disorder, whereas the width is not, since the 
disorder controls relative weights of central and periph- 
eral regions of the broaden ZBCP. For a weak broadening 



and high experimental resolution impurity bands can be 
identified, at least for low ns and u, as low-energy humps 
near the zero-energy peak. When Ug goes up, the bands 
become more close to the zero energy. Then more resolu- 
tion and less broadening are needed for identifying them 
as separated from the zero-energy peak. For concentra- 
tions ns > n-c impurity bands with positive and nega- 
tive energies merge into one impurity band, centered at 
zero energy. Sources, not related to the surface disorder, 
transform an abrupt change of the regimes of broadening 
of the ZBCP, into a gradual crossover from one regime to 
another. While the height of the ZBCP can be dominated 
by the surface disorder for any ns, its width is always in- 
dependent of ns below the crossover region. Above the 
crossover the surface disorder can strongly influence the 
broadening of the ZBCP. It is worth noting, that special 
experimental study of effects of surface disorder on the 
shape of the ZBCP has been carried out in Ref.l^^ where 
various degrees of the surface quality have been prepared 
by changing ion irradiation of YBa2Cu307_5/Pb junc- 
tions. It was found in a certain range of irradiation doses, 
that the ZBCP width is almost unchanged as a function 
of surface disorder, while the height of the peak effec- 
tively decreases with increasing the irradiation effects. 
This kind of behavior is in agreement with theoretical 
results of the present paper. 

We thank P. Hirschfeld and N. B. Kopnin for useful 
discussions. This work was supported, in part, by the 
Russian Foundation for Basic Research under Grant No. 
02-02-16643 (M.S.K. and Yu.S.B.) and by the Swedish 
Rcseach Council (VR)(M.F.). 



1 S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 63, 1641 
(2000). 

^ T. Lofwander, V. S. Shumeiko and G. Wendin, Super- 

cond. Sci. Technol. 14, R53 (2001). 
^ A. V. Balatsky, M. I. Salkola, and A. Rosengren, 

Phys. Rev. B 51, 15547 (1995). 

* M. I. Salkola, A. V. Balatsky, and D. J. Scalapino, 
Phys. Rev. Lett. 77, 1841 (1996). 

^ Yu. S. Barash and A. A. Svidzinsky, JETP 84, 619 (1997) 

[Zh. Eksp. Teor. Fiz. Ill, 1120 (1997)]. 
® Yu. S. Barash, Phys. Rev. B 61, 678 (2000). 

N. M. Chtchelkacthev and Yu. Nazarov, Phys. Rev. Lett. 

90, art. 226806 (2003). 

* J. Bardeen, R. Kiimmel, A. E. Jacobs, and L. Tewordt, 
Phys. Rev. 187, 556 (1969). 

® F. J. Culetto, G. Kieselmann and D. Rainer in Proceedings 
of the 17th International Conference on Low Temperature 
Physics, edited by U. Eckern, A. Schmid, W. Weber and 
H. Wiihl (North Holland, Amsterdam, 1984), p. 1027. 

1° Yu. N. Ovchinnikov, Sov. Phys. JETP 29, 853 (1969) [Zh. 
Eksp. Teor. Fiz. 56, 1590 (1969)]. 

" N. Schopohl and K. Maki, Phys. Rev. B 52, 490 (1995); 
N. Schopohl, cond-mat/9804064 (unpublished). 
M. Eschrig, Phys. Rev. B 61, 9061 (2000). 
E. M. Baskin, L. N. Magarill, and M. V. Entin, 
Zh. Eksp. Teor. Fiz. 75, 723 (1978) [Sov. Phys. JETP 48, 



365 (1978)]. 

A. M. Dyugaev, P. D. Grigoriev, and Yu. N. Ovchinnikov, 
Pis'ma Zh. Eksp. Teor. Fiz. 78, 180 (2003) [JETP Letters 
78, 148 (2003)]. 

I. S. Burmistrov and M. A. Skvortsov, Pis'ma 
Zh. Eksp. Teor. Fiz. 78, 188 (2003) [JETP Letters 78, 
156 (2003)]. 

A. Poenicke, Yu. S. Barash, C. Bruder, V. Istyukov, Phys. 
Rev. B 59, 7102 (1999). 

The presence of well pronounced zero-energy peak in the 
surface LDOS does not inevitably mean the large height of 
the ZBCP. For thick junctions only active tunneling tra- 
jectories, focused within a narrow cone around the surface 
normal, contribute to the conductance. This can strongly 
reduce the ZBCP (as well as the impurity peaks), with 
respect to the background level. 

Analogous condition can be obtained for scatterers, which 
are homogeneously spread out as bulk impurities in a d- 
wave superconductor up to (110) smooth surface, assum- 
ing that the zero-energy peak in the LDOS is sufficiently 
high. In terms of results of Ref. we find in this case, 
that the condition for the Born approximation to apply to 
low-energy quasiparticle scattering by impurities near the 
surface takes the form n 2> A''/Ao. 

Yu. Barash, H. Burkhardt, A. Svidzinsky, Phys. Rev. B 
55, 15282 (1997). 



9 



Since the TDL considers only the Born scatterers, the con- 
dition ns ^ ria concerns, in particular, the applicability of 
the results for effects of the TDL on the low-energy spectra 
of d-wave superconductors. 

^° A. Poenicke, M. Fogelstrom, J. A. Sauls, Physica B 284- 
288, 589 (2000). 

21 P. J. Hirschfeld, N. Goldenfeld, Phys. Rev. B 48, 4219 
(1993). 



N. B. Kopnin, Phys. Rev. B 65, 132503 (2002). 
M. B. Walker and P. Pairor, Phys. Rev. B 59, 1421 (1999). 
T. Lofwander, M. Fogelstrom, and J. A. Sauls, Phys. Rev. 
B 68, 054504 (2003). 

M. Aprili, M. Covington, E. Paraoanu, B. Niedermeier, 
and L. H. Greene, Phys. Rev. B 57, R8139 (1998). 



